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The formalism to describe electron scattering reactions on two-neutron halo nuclei is developed. 
The halo nucleus is described as a three-body system (core+n+n) , and the wave function is obtained 
by solving the Faddeev equations in coordinate space. We discuss elastic and quasielastic scattering 
using the impulse approximation to describe the reaction mechanism. We apply the method to 
investigate the case of electron scattering on ^He. Spectral functions, response functions, and 
differential cross sections are calculated for both neutron knockout and a knockout by the electron. 

PACS numbers: 25.30. Fj, 25.60.-t, 25. 70. Be, 21.45. +v 
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I. INTRODUCTION 

In the middle of the 80 's it was experimentally found 
that some light nuclei close to the neutron drip line have 
an spatial extension much larger than expected according 
to its mass number [|]-Q . This is the case for instance of 
^^Li, that has a root mean square radius similar to a nu- 
cleus with a mass number three times larger. Very soon 
it was suggested that these nuclei could be understood as 
a core surrounded by one or more neutrons, which extend 
to several times the nuclear radius ||5|. This picture was 
supported by subsequent studies that proved the valid- 
ity of the few-body models as an appropriate method to 
describe the most general properties of halo nuclei [p|-p^ . 

The term halo nuclei has then been coined to de- 
scribe very weakly bound and spatially extended systems, 
where several (usually one or two) neutrons have a high 
probability of being at distances larger than the typi- 
cal nuclear radius. For a general overview of their basic 
properties see for instance [|ll|. Among halo nuclei a lot 
of attention has been paid to two-neutron halo systems 
from which ^He (a -I- n -I- n) and -'^-'^Li (^Li +n + n) are 
their maximum exponents |12| . They are samples of the 
so called Borromean nuclei, that have the property of be- 
ing bound while all the two-body subsystems made by 
two of the three constituents are unbound. 

The ability to produce secondary beams of halo nuclei 
opened the possibility of investigating their structure by 
the measurement of the momentum distributions of the 
fragments coming out after a collision with stable nuclei 
[0-|l9|. The simplest approach to the understanding of 
Borromean nuclei fragmentation reactions was made by 
means of the sudden approximation P0|-p2| , that assumes 
that one of the three particles in the projectile is instan- 
taneously removed by the target, while the other two 
particles remain undisturbed. Clearly this can only be 
justified for reaction times much shorter than the char- 
acteristic time for the motion of the three particles in 
the system. Since the system is weakly bound this re- 
quirement is well fulfilled for a high energy beam. This 
model was proved to be valid to describe the core mo- 
mentum distributions p2| , but it failed in the attempt 
of reproducing experimental momentum distributions for 



neutrons. Several authors suggested that the final inter- 
action between the two non-disturbed particles played an 
essential ro le, esp ecially when low-lying resonances are 
present |l^j2l|,^. In corpor ation into the model of the 



final state interaction 1 24 25[| is in fact necessary in order 
to obtain a good agreement between theory and experi- 
ment also for neutron momentum distributions. Indeed 
further refinements P6| , P7| in the description of the reac- 
tion process have been proved to be crucial for a better 
interpretation of the experimental data. 

Although a lot of information has been extracted from 
the halo nucleus-nucleus fragmentation reactions, it may 
pay to look for cleaner ways of investigating the struc- 
ture of halo nuclei. As in any nucleus-nucleus collision 
the effects of the strong interaction involved in the re- 
action mechanism are mixed with those determining the 
nuclear structure. It is very well known that the way to 
avoid this problem is to substitute the hadron probe by 
electrons. In fact electron scattering is usually considered 
as the most powerful tool to investigate nuclear structure 
p8| , p9[ . The electromagnetic coupling constant is small 
enough such that the Born approximation can be used, 
the electron-nucleus interaction is well described through 
QED, and it is possible to vary independently the mo- 
mentum transfer and the energy transfer to the nucleus. 
The problem when dealing with halo nuclei is that they 
are far from the stability region, and the nuclear target 
can not be at rest as in conventional nuclear structure 
studies with electron beams. However it is in principle 
possible to perform electron scattering experiments by 
means of the collision between a secondary beam of un- 
stable nuclei and an electron beam. This kind of experi- 
ments are part of the MUSES (MultiUSse Experimental 
Storage rings) project in RIKEN, and they are projected 
for the beginning of the coming century [ pO[ . 

In parallel with the experimental projects it is then 
clear that theoretical studies of the process should be 
developed. From the experimental point of view the 
simplest reaction would be the elastic scattering. This 
would permit to investigate the charge density of the 
halo nucleus. This is specially interesting when the pro- 
ton dripline is approached because it directly gives the 
extension of the possible proton halo (*B is a good can- 



didate for this). For neutron halos the charge radius is 
expected to be the one of the core. Elastic electron scat- 
tering on ^^Li or ^He will be important to confirm that 
this picture of neutron halo nuclei is correct. As one 
slightly increases the energy transfer to the nucleus the 
halo will break up, and coincidence (e, e'x) measurements 
can be made, with x either a halo neutron or the stable 
core. 

The main goal of this paper is to contribute to the 
first steps in the investigation of electron scattering reac- 
tions on two-neutron halo nuclei. To illustrate the model 
we will consider the case of electron-^He scattering. Al- 



though estimates of the elastic scattering reaction will be 
given, most of the work will be devoted to (e, e'n) and 
(e, e'a) reactions at quasielastic kinematics. 

The paper is organized as follows: A brief descrip- 
tion of kinematics and of elastic scattering reactions is 
given in sections II and III, respectively. The differential 
cross section for exclusive (e, e'x) processes is given in 
section IV. Sections V and VI are devoted to the theo- 
retical description of the different ingredients involved in 
the electron-^He scattering process and to the presenta- 
tion of the results, respectively. Concluding remarks are 
given in section VII. 
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FIG. 1. (a) Scheme of a general electron-nucleus reaction in the Born approximation, (b) Axes system chosen to describe 
the process. 



II. KINEMATICS AND GENERAL 
CONSIDERATIONS 

In fig. la we show the scheme of a general electron scat- 
tering reaction. An electron with energy and momen- 
tum (EotPq) hits a target with energy and momentum 
{EfnPfi). The energy and momentum transfer to the nu- 
cleus is denoted by (w, q), and the energy and momentum 
of the scattered electron and the final nuclear system are 
given by (E'^^Pq) and {E'f^,p'^), respectively. 

In fig. lb we show our axis system. The z-axis is chosen 
along the momentum transfer q. The y-axis is defined 
along Pq X Pq, and it is therefore perpendicular to the 
scattering plane defined by pg and Pg . The x-axis is de- 
fined by y X 2:, where z and ij are unit vectors along the 
z-axis and y-axis, respectively. The angle between pg 
and Pq is the scattering angle 9e ■ 

Energy and momentum conservation in the process de- 
termine that 



Eo- 
Po 



Eh — Eq 



Ei 



- Ph =Po+Ph 



(1) 



and the energy and momentum transfer are related to 
the electron and nuclear energy and momentum by 



Lo — Eq — Eq — E,^ — Eh 
Q = Po-Po= P'h - Ph 



(2) 



Working in the frame of the nucleus {Eh — Mh and 
p,j = 0) these expressions can be rewritten as 

uj + Mh = K; q = p'h (3) 

and the invariant hadronic mass {W) takes the form 

W^^E';^-p'i = {u + Mhf-q^ (4) 



From eq.(|2|) we obtain that 



f = -ql = q^-uj^ = AEQE'Qs\n^^ 



(5) 



where ultrarelativistic electrons have been assumed (the 
electron mass is neglected). 

In the case of elastic scattering the invariant hadronic 
mass coincides with the mass of the target {W = Mh), 
and eq.(Q) leads to 

uj = {Ml + g2)i/2 _ Mh = E',^ - Mh = n (6) 

that is the kinetic energy of the final nucleus: 

/02 



L^ = Th = 



2Mh 



(7) 



Assuming that the two-neutron halo nucleus has no 
excited bound states (as it happens in ^He) an increase 



of the energy transfer to the nucleus will break the halo 
system into its three constituents. The energy and mo- 
mentum of the final hadronic system can then be written 
as 



Bh 



P? 



P2 






P'h = Pi + P2 + P3 



(8) 



where {Ei,p'^) {i — 1,2,3) are the energy and momen- 
tum of the three fragments, eqs. (|l|) to (Q) are then 
valid in this region after substitution on them of eq. <M). 
Processes involving excitations and one nucleon knock- 
out from the core may take place at higher energies, and 
will not be discussed here. 
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FIG. 2. Scheme of the reaction and definition of the co- 
ordinates used to describe the exclusive (e, e'x) process on a 
two-neutron halo nucleus. Only one of the constituents of the 
halo nucleus is considered to interact with the electron. 

These breakup processes will be studied assuming that 
the whole energy and momentum transfer are absorbed 
by one of the constituents (participant), which is knocked 
out from the halo system. The other two constituents are 
mere spectators. The large extension of the nucleus sup- 
ports this picture, since the probability for simultaneous 
interaction with the electron of two of the nuclear con- 
stituents is small. The scheme of the reaction can be seen 
in fig. 2. If we consider constituent 3 as the participant 
one, we then have 



E' 



P3 = P3 + 9 

Eh — E^ = Ei+ E2 

Ph-P3= P'l + P'l 



(9) 



where (£^3,^3) are the energy and momentum of the par- 
ticipant constituent inside the halo nucleus. 

In the target frame {Eh = ^h and p,j = 0) we have 
that pg = —(p'l +P2); and 

w -I- Mh ^E'^^m^+ 7712 + ma + Ti + Ta + T3. (10) 

Assuming non relativistic kinetic energies for the three 
fragments we have 



2m\ 2m2 27713 



EL 



vl 



(P3 



2(mi +7772) 



27773 



(11) 



where Bh is the halo binding energy, mi and Ti are the 
masses and final kinetic energies of the three constituents 
of the halo nucleus, and E'^ is the kinetic energy of the 
system made by particles 1 and 2 referred to its own 
center of mass. 

From eq. (^l|) one can see that for a fixed value of the 
momentum transfer q the variation in the energy trans- 
fer u) is connected to the different values of the internal 
momentum pg of the participant particle. Actually, for a 
fixed value of the momentum transfer, electron-nucleus 
scattering cross sections show a broad peak known as the 
quasielastic peak whose center is placed at pg = 0. Ne- 
glecting the binding energy Bh and the internal energy 
E'^ we obtain from eq.(O) the following value for the 
energy transfer in the center of the quasielastic peak 



(12) 



27773 



that corresponds to the kinetic energy transferred to a 
particle of mass 777g at rest. Reactions taking place in 
a quasielastic peak region are interpreted as processes 
where a single constituent in the nucleus absorbs the 
whole energy and momentum transfer, being ejected from 
the nucleus. 

When the participant particle is the core, the center of 
the quasielastic peak is at smaller energy than when the 
participant particle is one of the halo neutrons. In the 
^He case Wq.p. is approximately four times smaller for a 
participant ^He than for a participant neutron. 

In this paper the exclusive (e, e' x) reactions will be 
investigated in the so called perpendicular kinematics^ 
where the energy and momentum transfer are maintained 
fixed. In this kinematics, for a given electron beam with 
energy and momentum (_Eo,Po), the scattering angle 0^ 
and the energy of the outgoing electron Eq are fixed dur- 
ing the experiment. From (g) and (|^) it is clear that 6^ 
and Eq can be chosen in such a way that uj and q take 
the desired values. In particular, for a fixed value of q 
the energy transfer will be taken as shown in (|l2|) . We 
are then working in the center of a quasielastic peak. 

The energy transfer in the center of the quasielastic 
peak obviously increases with the momentum transfer q. 
If this momentum is such that uj is larger than the sepa- 
ration energy of the nucleons in core, then a neutron de- 
tected in coincidence with the electron could come from 
the core instead of coming from the halo. Therefore if 
Lu is maintained below the separation energy of the nu- 
cleons in the core we guarantee that the halo nucleus is 
broken in only its three constituents. In particular, for 
^He, since the neutron or proton separation energy in He 
is around 20 MeV a momentum transfer smaller than 200 



MeV/c would satisfy this condition. At the same time a 
value of uj below 20 MeV is also unable to excite the a 
particle to its first excited state. 

From eq.(^ and assuming that E^ = p\ 
write 



(p^ 



2p3gcos( 



2\l/2 



{pI + '^ 



2\l/2 



(13) 



If we now take q = 200 MeV/c and let the energy transfer 
take the value q^ /2ma {uj ^ 5 MeV), it is then simple to 
see from eq.(n3) that this energy can be absorbed by a 
participant a particle {m^ = ma) with rather low inter- 
nal momentum (smaller than a few MeV/c). If for the 
same value of uj we consider that the participant parti- 
cle is one of the halo neutrons (ma = m„) the required 
internal momentum of the neutron has to be very large 
(typically of several hundreds of MeV/c). However the 
probability of finding a halo neutron in ^He with such a 
large momentum is basically zero pl| , p6| . In the same 
way, from eq.(|l3|) one sees that if uj = q^/2mn (uj ~ 20 
MeV for q = 200 MeV/c) a participant neutron with 
small momentum can absorb such an energy, but it can 
not be absorbed by a participant a particle. 

Therefore, assuming that the reaction picture shown in 
fig. 2 is valid (a single constituent acts as participant while 
the other two are just spectators), we can then conclude 
that an energy transfer in the center of the a quasielastic 
peak (uj = q^ /2nia) selects the a as participant particle, 
and a halo neutron if uj is in the center of the neutron 
quasielastic peak {uj = q^ /2mn). 



III. ELASTIC SCATTERING 

The general expression for the elastic electron scatter- 
ing cross section can be found for instance in |Q. In 
particular for a nucleus with spin zero, as ^He, the elas- 
tic cross section in the nucleus frame takes the form 



da Z'^aM 



dCle 



Ire 



\Fchiq)\' 



(14) 



where the ultrarelativistic limit for the electrons has been 
assumed, fte defines the direction of the outgoing elec- 
tron, and 

a^cos^de/2) 2EoS^n{(h/2) ,^ _, 



iE^ Sin'' (6^/2)' 



Mh 



are the Mott cross section and the recoil factor. 
The charge form factor Fch{q) is given by 

, , 47r /■ sin((7r) , . 2 . 

Fch{q) = -^ / ^-^p{r)r'dr 

Z J qr 



(16) 



where p(r) is the charge density. 

Assuming that the picture of a core plus several halo 
neutrons is valid to describe neutron halo nuclei, the 



charge density should be practically identical to the one 
of the core. Therefore, since the neutron contribution to 
the charge form factor will be negligible at moderate q- 
values, the charge form factor should be the one obtained 
in an elastic electron-core scattering process. 
1 
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FIG. 3. Charge form factor after elastic electron scatter- 
ing from He assuming that the charge density corresponds 
to the a charge density (solid line) , and assuming that charge 
is spread out over the whole ^He nucleus (dashed line). The 
inset shows the charge density used in both cases. 

To simplify the analysis we take the experimental 
charge density of '^He that can be found in table V of 
p2| . The root mean square radius of the charge distri- 
bution is 1.68 fm, and the charge form factor is shown 
in fig. 3 by the solid line. This should then be the charge 
form factor obtained after elastic electron scattering on 
^He provided that the three-body picture is correct. 

The other limit would be to consider that the charge is 
not contained only in the core but distributed like mat- 
ter over the whole nucleus. In that case the charge form 
factor would be the one shown in fig. 3 by the dashed line. 
The form factor is now narrower due to the larger value 
of the charge radius. This curve has been obtained by 
modifying the range of the gaussians describing the "^He 
charge density in [p2| in order to get a root mean square 
radius of the charge distribution of 2.5 fm, similar to the 
size of ^He. 

From fig. 3 we can conclude that elastic electron scat- 
tering from ^He, or in general from neutron halo nuclei, 
is an excellent tool to determine to what extent the few- 
body description is valid. The inset in fig. 3 shows the 
charge density for two cases mentioned above. 

Center of mass corrections, also contained in the recoil 
factor of eq. ( p^ , may slightly modify the form factor and 
the charge r.m.s. radius of ^He relative to that of ^He. 
These effects could be accounted for using the three-body 
wave functions employed in the next sections. However, 
to give a first estimate of the size of the differential cross 
sections the two extreme pictures presented in fig. 3 serve 
the purpose. For instance, for Eq « 200 MeV and 6e « 30 
degrees, the differential cross section is of around 100 pb. 



IV. QUASIELASTIC SCATTERING 

In this section we enter in the study of the coincidence 
(e, e'x) processes, with x either a halo neutron or the 
core. The scheme of the reaction was shown in fig. 2. A 
three-body system with energy and momentum {Eh,Ph) 
interacts with an electron with energy and momentum 
(Eo,Pq). The whole energy and momentum transfer 
(w, q) is absorbed by one of the three constituents of the 
nucleus (constituent 3) that is ejected out from the three- 
body halo system. The energies and momenta in the final 
state are denoted by {E[,p[) where i — Q for the electron 
and i = 1, 2, 3 for the three halo constituents. 

Together with the impulse approximation (the virtual 
photon is absorbed by a single constituent) we also ne- 
glect the Coulomb distortion of the ingoing and outgoing 
electrons. They are then described as plane waves. At 
the same time the distortion due to the final strong inter- 
action between particle 3 (participant) and the two-body 
subsystem made by particles 1 and 2 (spectators) is also 
neglected, and the knocked out constituent is also de- 
scribed as a plane wave. On the contrary, the interaction 
between particles 1 and 2 is included in the calculation. 
This interaction is what in refs. |2jj2^ is referred to as 
final state interaction, and plays an essential role in the 
description of the halo nucleus fragmentation reactions. 
Our picture is then similar to that described for instance 
in ||,|1 for (e, e' N) reactions in PWIA. The only dif- 
ferencefbut important) is that our residual nucleus is an 
unbound two-body system, while they consider that the 
residual nucleus is in a well defined bound state. 

In appendix A it is shown that working in the frame of 
the three-body halo nucleus (p/j = Pi +P2 +P3 = 0) the 
transition matrix element of the reaction has the form 



(^) 



Sf^ = T73 ^(2^)"^''^(^^ - Ef)5\P^ - Pf) 



(17) 



i(Po> cro)7pw(Po> c^o) Yl K'aM'P-i^^^'^^^<yi2,s^<yAP'x^P'y) 



E'2 



E'^ are 



where Ei ^ Eq + Eh and Ef = E'q + E[ 

the initial and final total energies, and analogously for 

the initial and final total momenta Pi and Pf. The Ja- 



cobi coordinates x and y (see fig. 2) are defined in (Al) 
and their conjugated momenta p^ and p are given in 
(A2). The momenta in the final state are denoted with 
primes. The spinor u(p, a) describes a free electron with 
momentum p and third spin component cr. J'^, (QjPs) 
is the matrix element in momentum space of the current 
operator connecting st ates o f the parti cle 3 with spin pro- 
jections (73 and (T3 (eq.(A16)). In (eq.(Al7)) we have also 
defined M"^^^ , (p',p'), which is the overlap between 

the initial three-body halo wave function ($-^^^(a;,y)) 
and the continuum wave function of the final two-body 
subsystem {w''^^''^^'^{x,p'^)). When the interaction be- 



tween particles 1 and 2 is neglected M, 



JM 

S120'12,S3(t', 



iPx.P'v) 



becomes the Fourier transform (normalized to 1) of the 
three-body halo wave function. S12 and ai2 are the to- 
tal spin and its third component of the two-body system 
made by particles 1 and 2. 

We can now compute the differential cross section 
as the transition probability per unit volume and unit 
time {\Sfi\'^/VT) divided by the fiux of incident parti- 
cles {po/VEq in the halo nucleus frame), divided by the 
number of target particles per unit volume (1/1^) and 
multiplied by the number of final states that is given by 



y dp'o ^^ dp[ ^^ dp'^ ^^ dp'^ 



(27r)3 (27r)3 (27r)3 (27r)3 
y dp'o ^^ dp'3 ^^ dp'^ ^^ dp'y 



(27r)3 (27r)3 (27r)3 (27r)3 



(18) 



After averaging over initial states and summing over 
final states we obtain 

, ,/'!!,, , = {27T)^^SiE, - Ef)S{P, Pf) 
dp'odp'^dp'^dp'y q^ppo 

7; Y. I"(Po'^o)7m"(Po.^o)|' (19) 



2J 



ViSE 



Ma3 S120-12 



E'^"3...',('7-P3Xf.,.,s3.^(pL,p;: 



In this expression the integral over p' can be immedi- 
ately done by use of the del ta of momenta, getting that 
Py = 9 - P3 = -P3 (see eq.(|AT5|)). 



In the ultrarelativistic limit it is known that l35l 
^ Y l"(Po.c^o)7M'«(Po,cro)P = „p p, ??M^ (20) 



"^EqE'q 



where the leptonic tensor rj^i, is 

Vt^u = POf^p'o^ + Pol'P'q^, + -^qld,,,^ 

and g^i, is the diagonal metric (1,-1,— 1,-1). 
For a particle 3 with spin or 1/2 the matrix 



(21) 



M, 



*.JM 

Sl2<^12, ssas 



{p'x,P3)M;!Z,2,S3a'iP'x^P3) 



(22) 



is diagonal and has both diagonal terms equal. The dif- 
ferential cross section can then be written 



2J 






JM 
2(^12.53^3 



(P^,P3)| (23) 



Afo-3Si20-i2 



with a — e^/47r, and where we have defined the hadronic 
tensor 



^^^-2;^El<.^('?'P3) 



(24) 



that is the hadronic tensor associated to an electron- 
particle 3 scattering process. 

When particle 3 is one of the halo neutrons the 
electron- neutron cross section enters in (|2^), and has the 
standard form 



a^-{q,p,)^^^V,.R''' 



,(") , T/ -oH 



<JM [VLlZr + VtTZIP' + VltTI'-;;^ + VttTZ^tt) (25) 

where the kinematic factors V's and the response func- 
tions 7?.'s are given by 












Vtt 






7^(?^ = i?22 _ ^11 



(26) 



as it can be found for instance in |3qj37|| . 

If particle 3 is the a-particle the electron-a cross sec- 
tion reduces to 



2a^ E; 



% Eq 



(a) 



(27) 



Explicit expressions for 7?.'"^ and T?.*-"-* are given in the 
next section. 

Taking this into account the nine-differential cross sec- 
tion (123) takes now the form 



^^^^ = (2.)3J^.(i..-i.,).-(,,P3) 



E'o' 



2J 



— E 

+ 1 .. ^ 



M, 



,JM 

Sl2Cri2,S3CT3 



ip'.,P3)\ (28) 



MCT3S12CT12 
where the delta of energies imposes that 

E,-Ef^ Mh +Lo~ ip'^ + to2)1/2 _ E'^ 



2(mi + TO2) 



(toi -f 7712) = 0, 



Mh is the mass of the halo nucleus and 



EL 



2/X12 



(29) 



(30) 



is the internal energy of the final two-body system. 

Making use of the delta function one can integrate over 
P3, that leads to 



dEQtlflp' dilp' dp 

"^ E 



7 = (^TT) — jrec(J (q,P3) 



2 J 



M, 



Mh 



JM 



(p'x^Pa)] 



(31) 



A/0-3S12CT1 



where P3 = P3 — q, p'^ is obtained from (p9[), the angles 
flpi and Qp' define the directions of p/o and p'^, respec- 
tively, and the recoil factor free is given by 



Jrec — 



Lup'^j — E'^q cos 6p' 



MhP's 



(32) 



Finally, from ( |30| ) we have that dp'^ — p'^^i2dE'^dflp'^ , 
and then one gets the familiar expression 

d^a 



dE'^dVLp,^dVLp,^dE'^ 



(27r) — frecO- [q, p^jb [E^, p-^) [33) 



Mh 



where 



2J+1 ,, ^ 

A/(T3Si20-12 



df}„/ \M. 



JM 
si2cri2 ,530-3 



(pi,P3)|' (34) 



We have then arrived in eg. (|33|) to the usual factor- 
ization of the cross section for an (e, e'N) process in 
PWIA ||33 , p5| -p7| . The differential cross section is written 
as product of the cross section of an electron-particle 3 
scattering process and a function S{E'^,p^) called Spec- 
tral Function. The spectral function is interpreted as the 
probability for the electron to remove a particle from the 
nucleus with internal momentum p, leaving the residual 
system with internal energy E'^ [ p3[ . 

In stable nuclei the residual nucleus is usually consid- 
ered to be in a bound state, in such a way that its in- 
ternal energy can only take discrete values corresponding 
to the different excited states of the residual nucleus. If 
the missing energy (defined as the part of the initial en- 
ergy that is not transformed into kinetic energy) is kept 
within the appropriate limits it is possible to select a 
specific state for the residual nucleus. This amounts to 
fix the internal energy E'^ to the value corresponding to 
a given excited state. The spectral function then gives 
the momentum distribution of the particle such that af- 
ter removal leaves the residual nucleus in such an excited 
state. 

In our case the residual nucleus is not bound, but it is 
made by two particles flying together in the continuum. 
In other words, the internal energy E'^ = p'^ / (2^112) takes 
continuum values, and it is not possible to select a defi- 
nite state for the final two-body system made by particles 
1 and 2. Therefore integration over E'^ is required, and 



d^a 



dE'^dflp'^dflp'^ 



dEL 



d^a 



dE'^d^p'^d^p'^dE'^ 



(35) 



Note that the cross section <j^^ and the recoil factor 
depend on pg (pg = P3 — q), and pg is obtained from 
(29) that depends on E'^. Therefore, the integration in 
(35) involves the whole expression in (p3|), and the fac- 
torization disappears. 

Nevertheless, when the two-body system has reso- 
nances at low energy, E'^ can be neglected in eq.(p9|), and 
the value of pg, and therefore u'^^ and the recoil factor, 
are to a good approximation independent of E'^. Now 
the integral over E'^ involves only the spectral function, 
and we can then write 
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where 



that can be interpreted as the momentum distribution of 
the constituent particle 3. 

If we now introduce eq. (|25|) or (^7|) , we can finally write 
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for the cases of participant neutron and participant a, 
respectively. The TVs are the usual structure func- 
tions, and they are product of the momentum distribu- 
tion n(p3) and the corresponding response function VS^^ . 
The five-differential cross section (Q) is formally iden- 
tical to the expressions given for instance in p3| , p6| , p7| 
for the (e, e'K) cross section in PWIA assuming a resid- 
ual nucleus in a bound state. The response function for 
(e, e'a) will be specified later. 

Up to now we have always considered that the particle 
detected in coincidence with the electron is the partici- 
pant particle (particle 3). However the case where one 
of the spectator particles (particle 1 or 2) is detected in 
coincidence with the electron is also of interest. When 
one of the neutrons in the two-neutron halo nucleus is 
knocked out by the electron we actually have two neu- 
trons in the final state, one of them participant and the 
other one spectator. When one neutron is detected it is 
obviously not possible to know which one was detected, 
and therefore both cases should be considered. 



In eq. ( p^ ) we can use the delta of momenta to integrate 
over P3 instead of p^ (pg = q — p'y)- Following the same 
steps we then arrive to an expression analogous to 
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From eq.(A2) we obtain 
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(41) 



from which it is easy to see that dp'^dp'y — dp^dp'y 
(i-1,2). 

Using now the delta of energies we can integrate over 
p'y and then get 
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where i = 1,2 and p'y is obtained from (29) (using the 
relation p'y — q — p'3 and eq.(pl|)). The recoil factor free 
is in this case the value of the inverse of the derivative of 
( P9| ) with respect to p'y . 

From (E^) we obtain the five-differential cross section 
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(43) 



that is analogous to ( p5| ) and where the momentum p'^ 
(i = 1, 2) is now the momentum of one of the spectator 
particles in the final state. 

In principle the differential cross sections ( pSJ ) and (|42 



depend on five variables (-Eg and four angles). However 
electron scattering experiments are normally performed 
in what is known as in plane kinematics. This means 
that only particles in the scattering plane are detected. 
We then restrict ourselves to processes in the a;z-plane 
(see fig. lb), and the azimuthal angles ipp' , tpp' (in {^^), 
and (/3p' (in (HS)) are taken equal to zero. The number of 
variables is reduced now to three. Working in perpendic- 
ular kinematics the energy transfer to the nucleus is kept 
fixed, and for a given incident electron energy Eq one has 
Eq = Eq — oj. There are then only two variables left, 0p> 
and 6p' in 



51), and 9„i and 6„' in (43). Finally we also 



Pi 



know that pg = Pq — q, and multiplying this expression 
by Po we obtain in the ultrarelativistic limit 



cosOp'^ = 



En cos 9,, -E' 



(44) 



and therefore the differential cross sections in the in plane 
perpendicular kinematics are simply a function of the po- 
lar angle 0p' in ( |35| ) or 9p' in (p3|). 



As mentioned in the previous section, in case that we 
neglect E'^ in ( P9|) the expression (|3q ) is vahd, and the 
observables can be given as function either of 9p' or of the 
internal momentum of the participant particle p^. This 
is because 9p' and p^ are connected through the relation 
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where p'^ is obtained from 
glected. 



j) provided that E'^ is ne- 



SPECTRAL FUNCTION AND RESPONSE 
FUNCTIONS 



With this prescription the expressions for the structure 
functions in eg. (|26|) are given in appendix B. 

If the participant particle is the core (a-particle in our 
case) we also need the a^"" cross section. Since the a- 
particle has spin zero, only the component with /z = is 
non-zero in ( Al(: ) , and 



J'iq) = Piq) = j^ I dze^1-^p{z) 



(48) 



As a consequence only TZl in (26) is different from zero, 
and then 

CT^"=t7A/VL7^L=aM%|p(q + P3)l' (49) 



In order to obtain the coincidence differential cross 
sections given in (|3^) and (43), we need to specify the 
spectral function ( p4D of the halo nucleus and the ele- 
mentary electron-neutron responses, the electron-a re- 
sponses, and the cross sections. 

In particular, for the spectral function we need the 
probability function 
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(46) 



The calculation of this summation requires the knowl- 
edge of i) the intrinsic wave function of the three-body 
halo nucleus ^''^^(a;, y), ii) the continuum wave function 
y^s 12, <j 12 (^^^pi^-^ of the spectator syst em in the final state, 
and in) compute the overlap in eq.(Al7). When apply- 
ing the model to a particular case the calculation of the 
wave functions need in its turn the interactions between 
the three constituents in the halo system. In the case 
of ^He, the neutron-neutron and neutron-a interactions 
need to be specified. 



A. Electron— participant particle cross section 

When the participant particle is one of the neu- 
trons from the halo the electron-neutron cross section is 
needed. Electron scattering by a free nucleon is a well un- 
derstood reaction, and the only uncertainties come from 
the off-shell character of the nucleon interacting with the 
electron. For halo neutrons off-shell effects in the nu- 
cleon current can not be expected to be important, and 
we consider the CCl prescription 
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subject to the condition of current conservation (J^ 



7° , 



V?)- 



The charge density p{r) for the a-particle can be taken 
from |32] , where it is parameterized as a sum of gaussians, 
as discussed in section III. 



B. Intrinsic wave functions 

The three-body wave function of the halo nucleus is 
obtained by solving the Faddeev equations in coordinate 
space. The total wave function '^•^^ of the three-body 
system (with total spin J and projection M) is written 
as a sum of three components, each of them written as 
function of one of three possible sets of Jacobi coordi- 
nates [p9|. For each Jacobi set we construct the hy- 
perspherical coordinates (p, a, fix, ^y) defined in refs. 
||l2|,p9|. The volume element is given lay p^d^dp, where 
dfl ~ sin a cos^ adad^xdfly. For each hyperradius p 
\^JM jg expanded in a complete set of generalized angu- 
lar functions ^n {p,fli) 
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(50) 



where p~^/^ is related to the volume element, and the 
index i refers to the three sets of Jacobi coordinates. 

The angular functions satisfy the angular part of the 
three Faddeev equations: 



2m p^ " "jH^™ -^-^nj 
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where {i, j, k} is a cyclic permutation of {1, 2, 3}, to is an 
arbitrary normalization mass, Vjk is the interaction be- 
tween particles j and k, and A^ is the p-independent part 
of the kinetic energy operator. The analytic expressions 
for A^ and the kinetic energy operator can for instance 
be found in |39|. 

The radial expansion coefficients fn{p) are obtained 
from a coupled set of "radial" differential equations [p9J , 
i.e. 
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where the functions P and Q are defined as angular in- 
tegrals; 

3 

p„„,(p) ^ J2 ldm^:>{p,n)-^<i>[^}{p,n) , (53) 
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After obtaining the three-body wave function (BG) two 
of the components can be rotated to the third one, such 
that vj/-^*^ is written as function of a single set of Ja- 
cobi coordinates. In particular, to describe the reaction 
shown in fig. 2 it is convenient to write vj/-^^^ in terms of 
the Jacobi coordinates x and y shown in the same fig- 
ure (note that the Jacobi coordinates are usually defined 
with some mass factors ]39|] that for simplicity are not 



included in (Al)) 



The continuum wave function 



^^^{p'^,x) describ- 



ing the spectator two-body system in the final state is 
expanded in partial waves 

w''^^'''^^{p'^,x)^ (55) 

^ il2fl2mi2 

^ (^12TO/!i2;Sl2Cri2|jl2mi2)/^^Y«i2m,^^(l^p^) 
mi!i2=-^i2 

where the radial functions u''/^ ^ (p^, x) are obtained by 
solving the Schrodinger equation with the appropriate 
two-body potential p^l . £12 is the relative orbital angu- 
lar momentum between particles 1 and 2, that coupled 
to S12 gives the total angular momentum J12. When the 
final state interaction between the spectator particles is 
neglected the expansion ( |55| ) becomes the usual expan- 
sion of a plane wave in terms of spherical Bessel functions. 

C. Wave functions overlap and spectral function 

The cross sections derived in section IV are given in 
terms of the overlap between the initial two-neutron halo 
wave function and the continuum wave function of the 



spectator system in the final state (eq.( A17)). After writ- 
ing (^) in terms of the Jacobi coordinates defined in 
( Al) it can be shown that the analytic expression of this 
overlap is |23] 
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(56) 



k'y and a^ 



— arctan(fc^/fc' ) 



where k^ — k'^ 

K = V"VMl2Pi;7 k'y = y/m/pi2^3P3, Ml2,3 = (mi + 

1112)1113/ (mi + 7712 + 7773), and ttj is a normalization mass. 
£3 is the relative orbital angular momentum between the 
participant particle 3 and the center of mass of the spec- 
tator system 1-1-2, and it couples to S3 to give the angular 
momenta J3. The /-functions are given in eq.(16) of P5[ , 
and they are computed numerically. 

Inserting this into ( p4|) we obtain for the spectral func- 
tion the following analytic expression: 

S(E',,p,) = P^ y y n2'nLSL'S' (57) 

J12^12Sl2f3J3 LSL'S' 

{J J12 js] (j 312 J3I 
L £,2 4 U; ^12 ^3 
S S12 S3J [5" S12 S3J 

Note that the spectral function does not depend on the 
direction of Pg. 

D. Interactions 

To obtain the intrinsic wave functions (pOh and (K3) the 
interactions between the constituents of the halo nucleus 
have to be specified. 

As indicated in [Q the details of the radial shape of 
the neutron-neutron interaction are not very relevant for 
the ^He ground state wave function as long as the low en- 
ergy n—n scattering parameters are correct. We then use 
a simple potential that reproduces the experimental s- 
wave and p-wave scattering lengths and effective ranges. 
In particular we choose a potential with gaussian shape 
including a central, spin-orbit (£„„ • s„„), tensor (512), 
and spin-spin (si • S2) interactions. This potential was 
derived in BO], and it is given by 

Vnn(r) = 37.05 exp(-(r/l. 31)2) _ 7.38exp(-(r/1.84)2) 
s„„ + 7.16exp(-(r/2.43)2)5i2 
- 29.53 exp(-(r/1.84)2))si -82, 



-23.77 exp(-(r/1.45)2)£„ 
-f (49.40 exp(-(r/l. 31)2 



where .^„„ is the relative neutron-neutron orbital angular 
momentum and s„„ = Si + 82- 

For the neutron-a potential we take central and spin- 
orbit parts with gaussian shapes 



(58) 



V^%) = V^'\r) + V^,!,\r)lr, 



(59) 



where s„ is the spin of the neutron and £„„ is the rel- 
ative neutron-a orbital angular momentum. The pa- 
rameters of the gaussians are fitted to reproduce the 
phase shifts for s, p, and d waves up to 20 MeV. The 
gaussians used are [Q 48.2exp (— r^/2.33^) for the cen- 
tral s-wave potential, — 47.4exp (— r^/2.30^) for the cen- 
tral p-wave potential, and — 21.93 exp(—r^/2. 03^) for the 
central d-wave potential. For the spin-orbit interac- 
tion we take — 25.49 exp(—r^/l. 72'^) for all the waves. 
All the strengths are given in MeV and the ranges in 
fm. Together with the phase shifts this interaction re- 
produces the s-wave scattering length of —2.13 fm, and 
the two p-resonance energies and widths E{p^/2) — 0.77 
MeV, r(p3/2) = 0.64 MeV, ^(pi/a) = 1.97 MeV, and 
r(pi/2) = 5.22 MeV (see ^). 

Note that for the central s-wave interaction we use a 
repulsive potential. This is done to avoid the neutrons 
from the halo to occupy the s-states already occupied by 
the neutrons from the a-particle, and therefore forbidden 
by the Pauli principle. One could in principle have used 
an attractive potential and directly exclude in the calcu- 
lation the Pauli forbidden states. As far as the repulsive 
and attractive potentials have the same low-energy prop- 
erties they give almost indistinguishable wave functions 

It is well known as a general fact that two-body in- 
teractions accurately reproducing neutron-neutron and 
neutron-a low energy scattering data, systematically un- 
derbind the ''He system by roughly 500 keV Q. To 
alleviate this problem the range of all the interactions 
is normally increased by a few percent and the central 
p-wave potential is made a bit deeper |13M3]. However 
this procedure modifies the energy of the p-resonances, 
and many observables, as momentum distributions after 
fragmentation or invariant mass spectra, are very sensi- 
tive to the position of these resonances [£5|,^ . We have 
then preferred to maintain untouched the two-body in- 
teractions and introduce an effective three-body force in 
eq.(p2) for fine tuning. The idea is that the three-body 
force should account for the polarization of the particles 
beyond that described by the two-body interactions. All 
the three particles must be close to produce this addi- 
tional polarization, and therefore the three-body force 
must be of short range. In particular we use the attrac- 
tive gaussian -7.55exp (-p^/2.9^) Q. 

The two-body neutron-neutron and neutron-a inter- 
actions, together with the effective three-body force, pro- 
vide after solving the Faddeev equations a ^He binding 
energy of 0.95 MeV and a r.m.s. radius of 2.5 fm, in 



good agreement with the experimental data of 0.97±0.04 
MeV and 2.57 ± 0.10 fm, respectively. These results are 
obtained including s, p, and d waves in the calculation. 
Three terms in the expansion ( pOJ ) were included. How- 
ever 99.18% of the norm of ^'-^ was found to be given 
by the first term with n — 1, 0.75% is given by the second 
term with n = 2, and only 0.07% is given by the term 
with n = 3. Therefore only the first two terms are really 
needed in order to obtain an accurate enough ^He wave 
function. Exclusion of the n = 3-term does not modify 
the binding energy and r.m.s. radius. 

When the participant particle 3 is the a-particle the 
X coordinate connects the two spectator neutrons, and 
^12 is its relative orbital momentum. In the same way 
y connects the center of mass of the two neutrons and 
the a, and £3 is its relative momentum. After comput- 
ing the three-body halo wave function v]/-'^^ in terms of 
these {x,y) Jacobi coordinates it is seen that 87.5% of 
the norm is given by the s- waves, 9.9% by the p- waves, 
and 2.6% by the rf-waves. Therefore the relative mo- 
mentum between the two neutrons in ^He is basically 
an s-state, and the d-waves play a minor role. In the 
same way we can consider that one of the neutrons is the 
participant particle. In this case it is convenient to write 
^'-'*^ in terms of the {x, y) coordinates with x connecting 
the spectator neutron and the a-particle, and y connect- 
ing the participant neutron and the center of mass of the 
spectator system n + a. Doing this one sees that 10.8% of 
the norm is given by the s-waves, 88.5% by the p-waves, 
and 0.7% by the d-waves. Therefore the neutron and the 
a are mainly in a relative p-state. To be more precise 
82.3% of the norm is given by the p3/2-wave, and 6.2% 
by the Pi/2- All these data are summarized in table |. 

Due to the low d-wave content in the ^He wave func- 
tion only s and p waves will be considered in the compu- 
tations. The inclusion of the d-waves does not produce 
visible changes in the results. 

TABLE I. Contribution to the norm of the ^He wave func- 
tion from the s, p, and d waves with the interactions given 
in the text. The second column gives the contributions when 
the X coordinate connects the two neutrons and y goes from 
the center of mass of the two neutrons to the a-particle. The 
third column gives the contributions when the x coordinate 
connects one neutron and the a-particle and y goes from the 
center of mass of the n + a system to the second neutron. 
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FIG. 4. Spectral function (eq.(|34[)) after knockout of one of the halo neutrons in 
right side is the contour plot of the spectral function shown in the left part. 



'He by the electron. The figure in the 



VI. RESULTS 

In this section we give the results for the observables 
derived in section IV for an electron scattering process 
on ^He. The computations are done following the steps 
given in section V. 

All the observables can be computed in two different 
scenarios. First, assuming that the participant particle 
is one of the neutrons and the unbound ^He system sur- 
vives in the final state; and second, the a particle is the 
one knocked out by the electron and the two neutrons 
survive undisturbed in the final state. In both cases we 
consider that we are in the corresponding quasielastic 



peak, i.e. the energy transfer is a; = q^/2M where M is 
the neutron mass in the first case and the ''He mass in 
the second case. 



A. Spectral functions and momentum distributions 



The spectral function S{E'^^p^) is given in (34), while 
in (m) we show its analytic form obtained from the spe- 
cific expressions used for the wave functions. As already 
mentioned, this function gives the probability of remov- 
ing the participant particle with momentum p^ leaving 
the residual system with energy E'^. 
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FIG. 5. Same as fig. 4 for the case of a knockout by the electron. 



Let us start with the case of a participant neutron. 
The spectral function is shown in fig. 4, and reflects the 
momentum distribution of the internal nucleon and the 
energy spectrum of the unbound system ^He. In the left 
part of the figure we show the spectral function itself, 
while the right part shows the projection over the basis. 

The spectral function has a peak whose precise posi- 
tion is easily seen in the right part of the figure. The 
maximum value corresponds to a "'He energy of roughly 
E',^ = 0.75 MeV, that is the energy of the lowest P3/2 



resonance in ^He (0.77 MeV |42|). The width of this res- 
onance is relatively small (less than 0.7 MeV), and this 
makes the peak clearly pronounced. The next resonance 
in ^He (pi/2) has an energy of around 2 MeV, and it 
should also appear in the figure. However, its large width, 
more than 5 MeV, makes this resonance to be smeared 
out and hidden by the P3/2. Furthermore, as mentioned 
in subsection V.D, more than 82% of the norm of the ^He 
wave function is given by the Pz/2 wave, while the pi/2 
only gives around 6%. The presence of s-waves in the 
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^He wave function (^ 10%) makes the spectral function 
non zero along the p„ = edge. 

In fig. 5 we show the same as in fig. 4 but for the case 
where the a-particle is scattered away by the electron, 
and the two neutrons act as spectators. In this case the 
overlap (56) is computed writing the ^He wave function in 



terms of the Jacobi coordinates with x connecting both 
neutrons. From table || we know that both neutrons are 
preferably in a relative s-state, and this will produce a 
non zero spectral function along the Pa — axis. The 
fact that the spectral function is proportional to \/E'^ 
(see eq.(|5^)) will make the function equal to zero at the 
origin. This peaked spectral function close to the origin 
refiects the fact that the n — n interaction has a very 
large scattering length for the ^S'p-wave (more than 23 
fm) , giving rise to a virtual state at very low energy (less 
than 100 keV). 



participant neutron 



participant a 




100 

p„ (MeV/c) 



200 



100 

p„ (MeV/c) 



200 



FIG. 6. (a) Neutron momentum distribution (see eq.(p7|)) 
for the participant neutron case. The contributions to 
the total from the Si/2 wave (short-dashed), Pi/2 wave 
(dot-dashed), and P3/2 wave (long-dashed) are shown. The 
inset shows the same contributions multiplied by p^, where 
the weight of each wave can be better appreciated, (b) Same 
as (a) for the a momentum distribution in the participant a 
case. 

After integration of the spectral function over the inter- 
nal two-body energy E'^ we obtain the internal momen- 
tum distribution of the participant particle (eq.(|37|)). In 
fig. 6a we show the computed neutron momentum distri- 
bution after integration over E'^ of the spectral function 
in fig. 4 (solid line). The contributions to the total from 
the si/2 wave (short-dashed line), pi/2 wave (dot-dashed 
line), and P3/2 wave (long-dashed line) are also shown. 
It may look that the s-wave contribution is much larger 
than the 10% given in table I. This is because the whole 
contribution at the origin comes from the s waves. How- 
ever, when computing the normalization of the momen- 
tum distribution the phase volume contains a p^ factor, 
that reduces the weight of the s contribution. In the in- 
set of fig. 6a we show the neutron momentum distribution 
multiplied by p^. We clearly see that dominant contribu- 



tion comes from the P3 /2 wave, while the s contribution 
is small. 

In fig. 6b we show the same but for the case of partici- 
pant a. Again the solid line shows the total a momentum 
distribution, and the long-dashed and short-dashed cor- 
respond to the contributions from the s and p waves, re- 
spectively. In this case (table B) the dominant term is the 
relative neutron-neutron state with orbital angular mo- 
mentum zero. In fact the long-dashed curve and the solid 
one coincide almost in the whole range of a momenta. 
Even when the momentum distribution is multiplied by 
the p^ factor coming from the phase volume (inset in 
fig. 6b) it is seen that the s waves are the clearly domi- 
nant ones, and they are responsible for the pronounced 
peak observed in the spectral function in fig. 5. 



B. Cross sections and response functions 



To compute the cross sections ( p5| ) and (^) we work 
in the in plane perpendicular kinematics (see end of sec- 
tion IV) . The momentum transfer q is fixed at a value of 
q — 190 MeV/c, and the energy transfer corresponds to 
the center of the quasielastic peaks {uj = q^j2m3), with 
TO3 = TO„ or ?Ti3 — ma- The cross sections ( p5[) and (|43| ) 
are then functions of the polar angles of the particle de- 
tected in coincidence with the electron, i.e. 0p' and 9p' , 
respectively. 

i) Participant neutron. 

Let us first consider the case when u corresponds to the 
center of the neutron quasielastic peak. As discussed at 
the end of section II only a halo neutron can be the par- 
ticipant particle. If the particle detected in coincidence 
with the electron is the participant neutron ((e, e'n) pro- 
cess) the cross section is given by eq.(p5|). Assuming that 
the internal energy E'^ can be neglected in eq.(|2^) the 
cross section can then also be written as shown in (pq). 

In fig. 7 we show the differential cross section for this 
process. Two cases have been considered, one of them at 
a forward scattering angle {6e — 30) , and the second one 
at a backward scattering angle {de = 150). The electron 
beam energy Eq takes a value of 375 MeV in the first 
case, and 108 MeV in the second. In the figure the solid 
line is the calculation as given in eq. (^5|) . The dashed line 
is the calculation as in ( p8| ) where E'^ has been fully ne- 
glected in eq.(29), i.e. E'^ — 0. Finally the long dashed 
line is the calculation as in ( p8|) but taking E'^ = 0.77 
MeV in eq.(p9). This energy corresponds to the energy 
of the lowest p-resonance in ^He. 
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FIG. 7. Differential cross section for a (e, e'n) reaction 
for a participant neutron process. The participant neutron is 
assumed to be detected in coincidence with the electron. The 
calculations with scattering angles 9e = 30 and 6e — 150 are 
shown. The solid lines are calculations as given in (651). The 
short-dashed and long-dashed lines are calculations as given 
in (H) taking E'^ = and E'^ = 0.77 MeV in (||), respec- 
tively. The inset shows the structure functions appearing in 

As seen in the figure the three computed cross sec- 
tions are very similar, with only a small difference in the 
maximum, that is placed at a value of 9^/ of roughly 
10 degrees. When the approximation in ( p8| ) is used the 
chosen value for E'^ in eq.( p^) is not very relevant, and 
only a full calculation as in (|35|) reduces the maximum 



and the value at 9, 



P'n 



0. In all cases the cross section is 



already negligible for outgoing neutrons with polar angle 
larger than 75 degrees. 

The calculations with different scattering angles differ 
roughly in a global scale factor. This is because for a par- 
ticipant neutron only the magnetic nucleon form factor 
is important. Hence, in eq.(p8|) the main contribution 
is given by the transverse structure function, and the 
whole dependence on the scattering angle is contained 
in the global factor ctm^- From eqs.(|l5|) and ( [2^ ) one 
sees that (TmVt goes like 1/ sin^(6'e/2) and scattering an- 
gles close to zero give a much larger contribution than 
scattering angles close to w. This is seen in fig. 7, where 
the maximum value of the cross section is 8 times larger 
for 9e = 30 than for 6^ = 150. In the inset in fig. 7 we 
show the four response functions in eq. (p8[) . It is clear 
that, as mentioned above, the dominant contribution is 
by far given by the transverse part. Only Wlt gives an 
additional non negligible contribution. 

When the participant particle is one of the halo neu- 
trons, either the second halo neutron or the core can 
also be detected in coincidence with the electron. In this 
case the cross section is given by eq.(p^, where p'^ is 
the momentum either of the spectator neutron or of the 
spectator core (eq.(41)). The solid lines in fig. 8a give 
the differential cross section when the particle detected 
in coincidence with the electron is the second neutron. 
Therefore it also contributes to the e'n coincidence cross 



section. The two cases with 9e = 30 (thick lines) and 
9e = 150 (thin lines) are shown. As seen in the figure, 
for both scattering angles the maximum contribution to 
the cross section appears at an angle of around 70 de- 
grees, value for which the cross sections shown in fig. 7 are 
already small. Furthermore for small angles the contribu- 
tion to the e'n coincidence cross section coming from the 
spectator neutron is more than 10 times smaller than the 
contribution from the participant neutron (see fig. 7). In 
fig. 8b we show the total differential cross section (solid 
line) of the e'n coincidence reaction for 9e = 30. This 
curve is obtained by summing up the contribution from 
the participant neutron (short-dashed line in fig. 8b, or 
solid line for 9e = 30 in fig. 7) and the contribution from 
the spectator neutron (thick solid line in fig. 8a). From 
fig. 8b it is then clear that the behavior of the cross sec- 
tion is dominated by the participant neutron, while the 
spectator neutron contributes with a roughly constant 
cross section that creates a long tail that extends all the 
way up to 9p'^ — 180. When 9^ — 150 the difference is 
basically a global factor close to 8. We also show in fig. 8a 
the differential cross section for the e'a coincidence pro- 
cess when the a particle is a spectator. They are given 
by the dashed curves, and show a pronounced peak at 70 
degrees. 
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FIG. 8. (a) Differential cross sections for coincidence re- 
actions after a participant neutron process and assuming that 
the particle detected in coincidence with the electron is either 
the spectator neutron (solid line) or the spectator a (dashed 
line). Calculations with 9^ = 30 (thick lines) and 6^ — 150 
(think lines) are shown, (b) Total differential cross section 
after an electron-neutron coincidence reaction (solid line). 
The contributions from the participant neutron (short-dashed 
line) and the spectator neutron (long-dashed line) are shown. 
The dotted line is the differential cross section for a (e, e'a) 
process assuming that the a is the participant particle. 

a) Participant a. 

Let us finish with the case when tu corresponds to the 
center of the a quasielastic peak. The participant particle 
is then the a, m^ — nia, and for q = 190 MeV/c the en- 
ergy transfer is around 5 MeV. Neglecting E'^ in eq.(E9h, 
and taking M^ = nii + 77*2 + n^3, it can be seen that the 
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angle 9p' has to satisfy that sinOp' < (mi + 7712) /m^. 
This relation is fulfilled for any angle when the par- 
ticipant particle is one of the halo neutrons, but when 
TO3 > (ttii -\- 7712) it gives a limit to the value of 9p' . 
In particular, for a participant a particle one has that 
6pi ^ 30. This is also shown in fig. 8b, where the dotted 
curve shows the differential cross section for an exclu- 
sive (e,e'a) reaction in the center of the a quasielastic 
peak. The cross section is zero for angles larger than 
30, and the maximum value at the origin is around twice 
the one obtained in the e'n coincidence process. In the 
computation shown in the figure the scattering angle is 
30 degrees, for which Eq = 369 MeV {q = 190 MeV/c, 
uj — q^/2ma). As seen in eq.(|27|), only a longitudinal 
structure function enters in this case, and the whole de- 
pendence on the scattering angle is then contained in 
tJA/. Therefore different scattering angles only introduce 
a global factor determined by the ratio between the dif- 
ferent values of the Mott cross section. For instance, a 
scattering angle of 15 degrees would make the differential 
cross section four times larger than the one shown in the 
figure, while 6e = 150 would make the cross section more 
than 200 times smaller than the one in the figure. After 
a knockout by the electron the spectator neutrons could 
in principle be detected in coincidence with the outgoing 
electron. The computed cross section for this process is 
very small, and barely visible in the scale of fig. 8b. 



VII. SUMMARY AND FINAL REMARKS 

The fact that electron scattering is one of the most 
powerful and cleanest procedures in order to investigate 
nuclear structure lead us to consider these reactions to 
investigate one of the most intriguing objects recently 
discovered in nuclear physics: Halo nuclei. In partic- 
ular we have considered the case of two-neutron halo 
nuclei. The three-body wave function of the nucleus is 
obtained by solving the Faddeev equations in coordinate 
space by means of the adiabatic hyperspherical expan- 
sion. As soon as the two-body interactions between the 
constituents of the nucleus are specified the three-body 
structure can be computed to the needed accuracy. 

First we have made a short incursion into elastic elec- 
tron scattering reactions where we have seen the variation 
of the charge form factor according to two limiting values 
for the charge r.m.s. radius, the core size and the ^He 
size. The validity of the three-body picture requires that 
the charge form factor for ^He and "^He should be similar. 
In other words, if the present picture of ^He as a '^He core 
and two halo neutrons is correct, the elastic differential 
cross section for ^He should be practically identical to 
that of *He, while if the spatial charge distribution of 
^He would be as extended as that of the neutrons, the 
cross section would be at least four times smaller at typ- 
ical q values of 200 - 300 MeV/c. Thus elastic electron 
scattering provides a direct test of the halo picture. 



Then we make an extensive discussion of (e, e'x) pro- 
cesses, that are described in the impulse approximation, 
i.e. the virtual photon is absorbed by a single constituent 
of the halo nucleus. The quasielastic exclusive electron 
scattering formalism was developed years ago (see for in- 
stance [p3|), and in principle one could directly apply 
it to electron scattering on halo nuclei. However, the 
language used in the three-body description of the halo 
nucleus (inert core, Jacobi coordinates, ...) is not easily 
matched with the language commonly used in electron 
scattering, where the single-particle properties of nuclei 
play an important role. One of the main aims of this 
paper has been to obtain the differential cross section 
for electron scattering using the three-body wave func- 
tion of the halo nucleus as starting point. Comparison of 
the derived cross section with the well known expression 
for quasielastic exclusive electron scattering permits an 
easy identification of the observables of interest: Spec- 
tral function, response functions, and momentum distri- 
butions. 

After discussion of the kinematics of the reaction we 
have applied the method to electron scattering on ^He, 
for which the two-body interactions between its three 
constituents are well known. These interactions repro- 
duce the basic properties of the nucleus, separation en- 
ergy and root mean square radius, together with an im- 
portant amount of experimental data obtained after frag- 
mentation on stable targets. Uncertainties coming from 
a poor description of the halo nucleus are then mini- 
mized. We have investigated the different observables 
for the cases of a participant neutron and a participant 
a. 

We have first studied the spectral function, that carries 
the whole information about the structure of the two- 
body subsystem obtained after removal by the electron 
of one of the constituents of ^He. This function is in- 
terpreted as the probability of removing one of the con- 
stituents with a certain momentum leaving the residual 
system with a certain energy. Effects due to the presence 
of different partial waves can be observed, as well as in- 
dications about the resonance structure of the surviving 
unbound two-body subsystem. Of special interest is the 
case of neutron removal, since the spectral function di- 
rectly gives information about the unbound nucleus ^He. 

After integration of the spectral function over the en- 
ergy we obtain the momentum distribution of the partici- 
pant particle, that contains information about the weight 
of each partial wave in the wave function. Experimen- 
tal information about the momentum distribution can in 
principle be obtained by taking the ratio between the 
measured cross section and the electron-participant par- 
ticle cross section. However, one has to keep in mind 
that the factorization of the five-differential cross sec- 
tion is only approximated, and low lying resonances in 
the final unbound two-body system are required. For 
the same reason as before, low lying resonances in the 
two-body system are necessary in order to write the five- 
differential cross section in terms of the response func- 
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tions (that are simply the product of the momentum dis- 
tribution and the corresponding single-particle response 
function). This can always be done when the residual 
nuclear system is bound. We have seen that for ^He 
(that has a resonance at an energy of around 0.8 MeV) 
the approximated differential cross section does not differ 
very much from the exact one, and it is therefore possible 
to talk about response functions in this kind of reactions. 
Such a response functions have been obtained for the case 
of a participant neutron, where the transverse response 
clearly dominates, and for a participant a, where only 
the longitudinal one enters. 

In summary, as for any nucleus, electron scattering is 
probably the most accurate procedure in order to inves- 
tigate the structure of halo nuclei. Due to the up to now 
unattainable technical problems, experimental measure- 
ments are not available, and these reactions have been 
completely unexplored. Fortunately the first experiments 
of this kind are projected for the next few years [Q, and 
theoretical studies as the ones presented here are then 
needed. We have derived the cross section and investi- 
gated for the case of the two-neutron halo nucleus ^He 
how the observables are connected to the structure of 
the nucleus. We show results for differential cross sec- 
tions corresponding to elastic scattering as well as to 
quasielastic scattering on a neutron halo and on the ^He 
core, illustrating how the different kinematical regions 
probe the charge distribution, the spectral function and 
the momentum distributions in a transparent manner. 
Additional information could also be obtained exploit- 
ing spin polarization degrees of freedom, which are not 
discussed here given the incipient state of the subject. 



The coordinates x and y are the Jacobi coordinates 
except for the mass factors that are usually introduced 
and that for simplicity have been omitted here (see for 
instance [p5[). 

The conjugated momenta associated to the coordinates 
a;, y, and z are 



Px = , _ Pi \ P2 



mi + 7712 "nil + 7772 



P 



7773 



^ 7771 + 'T72 + TTil^ 



(Pl + P2) 



7771 + 1^2 

mi +7712 + 7773 



P3 (A2) 



Pz 



ms 7770 

-Po - _ , _ P3 



7770 +7773 mo + 7)73 



and the same with primes in the final state. The masses 
of the three halo constituents are denoted by 777^ (i — 
1,2,3). 

Working in the frame of the halo nucleus {pf^ — Pi + 
P2+ P3 = 0) one has 



Py=Pi+P2 = -Ps 



(A3) 



Following |44[| in the one photon exchange approxima- 
tion the transition matrix element for an electron scat- 
tering process is given by 

Sft = -ie / dtdroTpfit, ro)7pV'i(i, ro)A^'{t, Tq) (A4) 

where ^pij are the initial and final electron wave func- 
tions, that are given by 

Mt, ro) = -^n{po, ao)e-'^°*e'Po-^o (^5) 
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APPENDIX A: DERIVATION OF THE 
TRANSITION MATRIX 

In fig. 2 we show the scheme of the exclusive (e, e'x) 
reaction together with the notation for the different en- 
ergies and momenta involved. The axis system chosen to 
describe the process is given in fig. lb. If we denote by r^ 
(7 = 0, 1, 2, 3) the position of the electron and the three 
particles in the halo nucleus we then define the coordi- 
nates: 



ri -r2 

777iri -|- 77727-2 



mi + m2 
z = ro - r3 



-■^3 



(Al) 



TPfit, ro) = -^u{p',, a'o)e-^^'o'e^P'o-^o (A6) 



that are normalized to 1 in a volume V. (Tq and a'o are 
the third components of the spin for the ingoing and out- 
going electrons, respectively. 

The four-vector A^{t,ro) is the potential (scalar and 
vector) generated by the nucleus and seen by the electron 
at the instant t in the position ro. 



Substituting (|AJ) and (|A^) into ( |A4| ) one has 

Sr. = (A7) 

:^u{p'o,a'oh^uiPo, do) J dtdroe-'^'e'1-''»A'^it, ro) 

Making use of the Maxwell equations one has 

nA^(t,ro)-e(/|>|*) (A8) 

and the transition matrix element takes the form pq ] 

o2 



Sfi 



-76 



u{Po,<^ohMPo'<^o) 



^ fdtdroe~'^'e''l-''°{f\r\i) 

Hp J 



(A9) 
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where j^ is the nuclear current operator ( J'' = (p, J)). 

The initial hadronic state \i) is given by the three-body 
halo wave function that is written as 



with 



^•«^(a;^y)^e-^s,tg^p,.r. 



rh 



V 



miri + m2r2 + msri 

TOl + 7712 + 1^3 



(AlO) 



(All) 



and '^•^^ {x,y) is the intrinsic three-body halo wave 
function with spin J and third component M . 

The wave function of the two-body system made by 
particles 1 and 2 in the final state is written as 



1 
V 



w 



^i^''^i^(p;,a;)e-*(^;+^^)*e'(Pi+P^)-^i= (A12) 



where u!*^^''^^^(p^, a;) is the intrinsic continuum wave 
function, si2 is the coupling of the spins of both par- 
ticles and ai2 its third component, and ri2 — {miri + 

7772r"2)/(TOl +7^2). 

The fina l hadronic state |/) is made by the product 
of eq.(A12) and the plane wave describing the outgoing 
particle 3 



1 



-tEitm'^-rs 



(A13) 



where Xs3,(7' is the spin state of particle 3 in the final 
state. 

According to the scheme shown in fig. 2 particle 3 is the 
only constituent interacting with the electron. The rest 
of them are just spectators. The nuclear current operator 
should then be the one associated to the hadron labeled 
by 3. Assuming that this current operator depends only 
on the distan ce z b et ween particl e 3 an d the elec tron and 
substituting ( |A10D , ( |A12| ), and ( [aT^ ) into (^) a little 
algebra leads to 



nucleus, and then in the frame of the three-body halo 
system one has 

Pl+P2= Py = -P3 = Pi +P2=P'y (Al5) 

Therefore in this picture the center of mass momentum 
of the final two-body system made by particles 1 and 2 
reveals the internal momentum of the participant particle 
3. 

We denote now 

(A16) 

that is the matrix element in momentum space of the 
current operator connecting states of particle 3 with spin 
projections 173 and 0-3. 
We also define 



Si2CTi2,S3CT3 
, 1 



iP'x^P'y. 



' (27r)3 






Xss.a.,w^''-^^^-{x,p'J\^'^^'ix,y)) (A17) 



that is the overlap between the initial three-body halo 
wave function and the wave function of the final two- 
body system (the exponential is the two-body center of 
mass motion). When the interaction between particles 1 
and 2 is neglected the distorted function 7«*i2'°'i2(a;,p^) 



becomes a plane wave {e^-^^'^Xsi2.(Ti2) a-nd ( A17 ) is the 
Fourier transform (normalized to 1) of the three-body 
halo wave function. 



With (AI6) and (|A17|) we can finally write the transi- 
tion matrix element (|Al4) in the form 






u{Po,(^'ohp.uiPo,(^a) 



03 



J:..M^P3)MiX.,S3.M^P'y) (A18) 



^/^ = T73 72 (2'^)''5(^^ - Ef)5\P, ~ P/) (A14) 

u{p'Q,croh^,u{pQ,ao) 
J2 f dze^^-^{Xs3yjM^)\Xs3.'^3)i2nf 

T3 -^ 

{j^e'Py-yXs3^.sw'^^'^^H^,P'J\'^''H^,y)) 

where we have used that drf)dridr2dr^ = drodxdydz, 
and where the integrations over t and Tq give rise to 
the deltas forcing energy and momentum conservation. 
Ei^ Eo + Eh and Ef = E'q + E[+E'2 + E'^ are the initial 
and final total energies, and analogously for the initial 
and final total momenta Pi and Pf. 

In deriving ( [A14 ) we have also used that particle 3 ab- 
sorbs the whole energy and momentum transfer to the 



APPENDIX B: ELECTRON-NUCLEON 
OFF-SHELL STRUCTURE FUNCTIONS 

For the off-shell electron-nucleon cross section we have 
chosen what in refs. i39,pq| is called prescription CCl 
with current conservation. 

Let us call 



F,{ql) + F2{ql); 



(3- 



-F2{ql) 



2M 



nuclcon 



(Bl) 



The response functions defined as in eg. (pq) are then 
[in plane kinematics is assumed) 



7^L = 



2Ey,Er^ 



j[a^{2E:,E'^+ml-py^^)+ 



{Es + E'^f{P^{ml+py^^) + 2aPm^)\ (B2) 
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^T = in?r[«'(P?sin2 



£^3^3 






2v'i sin2 0p, (/32(^2 + p^p^^) + 2a/3m3)] (B3) 



7^ 



LT 



"^-p'^smOp-iEs+E',) 



E3E3 



(B4) 



TZtt — 
Finally 



1 



E3E2 



jp'i sin^ ^p, 



[a^ + 2p\ml +p^p'3^) + 4a/3m3] (B5) 



2 ^ G_E + tGm 



1 + T 



^2(<7^) = ^^f-^ (B6) 
^ 1 + r 



where r = — (?^/(4Afnucicon), and Ge and G^/ are param- 
eterized as in Bsl. 
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